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^ ■ Abstract 

O ■ We introduce the neural network approach to global fits of parton distribu- 

tion functions. First we review previous work on unbiased parametrizations 
of deep-inelastic structure functions with faithful estimation of their uncer- 
tainties, and then we summarize the current status of neural network parton 

O , distribution fits. 

■ We introduce the neural network approach to global fits of parton distribu- 

^ \ tion functions. First we review previous work on unbiased parametrizations 

J> ' of deep-inelastic structure functions with faithful estimation of their uncer- 

tainties, and then we summarize the current status of neural network parton 
^ , distribution fits. 

o\ : 

^ ' The requirements of precision physics at hadron colliders, as has been emphasized through this 

workshop, have recently led to a rapid improvement in the techniques for the determination of parton 
distribution functions (pdfs) of the nucleon. Specifically it is now mandatory to determine accurately the 
Q_i| uncertainty on these quantities, and the different collaborations performing global pdf analysis [1] have 

performed estimations of these uncertainties using a variety of techniques. The main difficulty is that one 
is trying to determine the uncertainty on a function, that is, a probability measure in a space of functions, 
and to extract it from a finite set of experimental data, a problem which is mathematically ill-posed. It 
is also known that the standard approach to global parton fits have several shortcomings: the bias intro- 
^ duced by choosing fixed functional forms to parametrize the parton distributions (also known as model 

dependence), the problems to assess faithfully the pdf uncertainties, the combination of inconsistent 
experiments, and the lack of general, process-independent error propagation techniques. Although the 
problem of quantifying the uncertainties in pdfs has seen a huge progress since its paramount importance 
was raised some years ago, until now no unambiguous conclusions have been obtained. 

In this contribution we present a novel strategy to address the problem of constructing unbiased 
parametrizations of parton distributions with a faithful estimation of their uncertainties, based on a com- 
bination of two techniques: Monte Carlo methods and neural networks. This strategy, introduced in [2,3], 
has been first implemented to address the marginally simpler problem of parametrizing deep-inelastic 
structure functions F(x,Q 2 ), which we briefly summarize now. In a first step we construct a Monte 
Carlo sampling of the experimental data (generating artificial data replicas), and then we train neural 
networks to each data replica, to construct a probability measure in the space of structure functions 
V [F(x, Q 2 )\. The probability measure constructed in this way contains all information from experi- 
mental data, including correlations, with the only assumption of smoothness. Expectation values and 
moments over this probability measure are then evaluated as averages over the trained network sample, 



F(x,Q 2 )})= / VFv\f(x,Q 2 )]t\f(x,Q 2 )} = — £ T (F^ k \x, Q 2 )) . (1) 



where T [F] is an arbitrary function of F(x, Q 2 ). 



The first step is the Monte Carlo sampling of experimental data, generating JV rep replicas of the 
original N^ at experimental data, 
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where r are gaussian random numbers with the same correlation as the respective uncertainties, and 
a stat , a sys ,aN are the statistical, systematic and normalization errors. The number of replicas N rep has 
to be large enough so that the replica sample reproduces central values, errors and correlations of the 
experimental data. 

The second step consists on training a neural network 1 on each of the data replicas. Neural net- 
works are specially suitable to parametrize parton distributions since they are unbiased, robust approxi- 
mants and interpolate between data points with the only assumption of smoothness. The neural network 
training consist on the minimization for each replica of the \ 2 defined with the inverse of the experimen- 
tal covariance matrix, 
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Our minimization strategy is based on Genetic Algorithms (introduced in [4]), which are specially suited 
for finding global minima in highly nonlinear minimization problems. 

The set of trained nets, once is validated through suitable statistical estimators, becomes the 
sought-for probability measure V [F(x, Q 2 )] in the space of structure functions. Now observables with 
errors and correlations can be computed from averages over this probability measure, using eq. Q. For 
example, the average and error of a structure function F(x,Q 2 ) at arbitrary (x,Q 2 ) can be computed as 

1 iVrep / 

F(x,Q 2 )) = — F^ k \x,Q 2 ), a(x,Q 2 ) = ^J(F(x,Q 2 ) 2 ) - (F(x,Q 2 )) 2 . (4) 

iVre P k =i 

A more detailed account of the application of the neural network approach to structure functions can be 
found in [3], which describes the most recent NNPDF parametrization of the proton structure function 2 . 

Hence this strategy can be used also to parametrize parton distributions, provided one now takes 
into account perturbative QCD evolution. Therefore we need to define a suitable evolution formalism. 
Since complex neural networks are not allowed, we must use the convolution theorem to evolve parton 
distributions in x— space using the inverse T(x) of the Mellin space evolution factor T(N), defined as 



q(N, Q 2 ) = q(N, Q 2 )T (n, a s (Q 2 ) , a s [q 2 



(5) 



The only subtlety is that the x-space evolution factor T{x) is a distribution, which must therefore be 
regulated at x = 1, yielding the final evolution equation, 



q(x, Q 2 ) = q(x, Q 2 ) £ dy T(y) + jf ^T(y) (q Q%j - yq(x, Q 



(6) 



where in the above equation q{x,Q 2 ) ) is parametrized using a neural network. At higher orders in per- 
turbation theory coefficient functions C{N) are introduced through a modified evolution factor, T(N) = 
T(N)C(N). We have benchmarked our evolution code with the Les Houches benchmark tables [5] at 

'For a more throughly description of neural network, see [2] and references therein 

2 The source code, driver program and graphical web interface for our structure function fits is available at 

|http : / / sophia .ecm.ub. es/ f 2neural| 




Fig. 1: Preliminary results for the NNPDF gjvs fit at Qq = 2 GeV 2 , and the prediction for _F^ vs (a:, Q 2 ) compared with the 
CTEQ and MRST results. 

NNLO up to an accuracy of 10 -5 . The evolution factor T(x) and its integral are computed and interpo- 
lated before the neural network training in order to have a faster fitting procedure. 

As a first application of our method, we have extracted the nonsinglet parton distribution qNs{x, Qq) = 
| (u + u — d — d) (x, Qq) from the nonsinglet structure function F^fx, Q 2 ) as measured by the NMC 
[6] and BCDMS [7, 8] collaborations. The preliminary results of a NLO fit with fully correlated uncer- 
tainties [9] can be seen in fig. ^ compared to other pdfs sets. Our preliminary results appear to point 
in the direction that the uncertainties at small x do not allow, provided the current experimental data, to 
determine if qNs{%, Q 2 ) grows at small x, as supported by different theoretical arguments as well as by 
other global parton fits. However, more work is still needed to confirm these results. Only additional 
nonsinglet structure function data at small x could settle in a definitive way this issue 3 . 

Summarizing, we have described a general technique to parametrize experimental data in an bias- 
free way with a faithful estimation of their uncertainties, which has been successfully applied to structure 
functions and that now is being implemented in the context of parton distribution. The next step will be 
to construct a full set of parton distributions from all available hard-scattering data using the strategy 
described in this contribution. 

References 

[1] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Nnlo global parton analysis, Phys. 
Lett. B531 (2002) 216-224, | |hep-ph/0201127| . 

D. Stump et at, Inclusive jet production, parton distributions, and the search for new physics, 
JHEP 10 (2003) 046, |hep-ph/ 03030 13| . 

S. Alekhin, Parton distributions from deep-inelastic scattering data, Phys. Rev. D68 (2003) 
014002, |hep-ph/0211096| . 

[2] S. Forte, L. Garrido, J. I. Latorre, and A. Piccione, Neural network parametrization of 
deep-inelastic structure functions, JHEP 05 (2002) 062, |hep-ph/020 42 32|. 

[3] NNPDF Collaboration, L. Del Debbio, S. Forte, J. I. Latorre, A. Piccione, and J. Rojo, Unbiased 
determination of the proton structure function f2(p) with faithful uncertainty estimation, JHEP 03 
(2005) 080, ||hep-ph/ 050 10 67| . 

3 Like the experimental low x deuteron structure function which would be measured in an hypothetical electron-deuteron 
run at HERA II, as it was pointed out during the workshop by M. Klein and C. Gwenlan 



[4] J. Rojo and J. I. Latorre, Neural network parametrization of spectral functions from hadronic tau 
decays and determination ofqcd vacuum condensates, JHEP 01 (2004) 055, 
|hep-ph/040 10471 . 

[5] W. Giele et ah, The qcd/sm working group: Summary report, hep- ph/02 04 316| 

[6] New Muon Collaboration, M. Arneodo et al, Measurement of the proton and deuteron structure 
functions, f2(p) andfl(d), and of the ratio sigma(l)/sigma(t), Nucl. Phys. B483 (1997) 3^-3, 
lhep-ph/9 6102 311 . 

[7] BCDMS Collaboration, A. C. Benvenuti et al.,A high statistics measurement of the proton 
structure functions f(2) (x, q**2) and rfrom deep inelastic muon scattering at high q**2, Phys. 
Lett. B223 (1989) 485. 

[8] BCDMS Collaboration, A. C. Benvenuti et al.,A high statistics measurement of the deuteron 

structure functions f2 (x, q**2) and rfrom deep inelastic muon scattering at high q**2, Phys. Lett. 
B237 (1990) 592. 

[9] NNPDF Collaboration, J. Rojo, L. Del Debbio, S. Forte, J. I. Latorre, and A. Piccione, The neural 
network approach to parton fitting, hep-ph/0 50 50 4 4| 



